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By direct numerical simulation of the time-dependent Gross-Pitaevskii equation using the split- 
step Fourier spectral method we study different aspects of the localization of a cigar-shaped inter- 
acting binary (two-component) Bose-Einstein condensate (BEC) in a one-dimensional bi-chromatic 
quasi-periodic optical-lattice potential, as used in a recent experiment on the localization of a BEC 
[Roati et al., Nature 453, 895 (2008)]. We consider two types of localized states: (i) when both local- 
ized components have a maximum of density at the origin x — 0, and (ii) when the first component 
has a maximum of density and the second a minimum of density at x = 0. In the non-interacting 
case the density profiles are symmetric around x = 0. We numerically study the breakdown of this 
symmetry due to inter-species and intra-species interaction acting on the two components. Where 
possible, we have compared the numerical results with a time-dependent variational analysis. We 
also demonstrate the stability of the localized symmetry-broken BEC states under small perturba- 
tion. 
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I. INTRODUCTION 

After the prediction of the localization of the electronic 
wave function in a disordered potential by Anderson fifty 
years ago [![, there have been many studies on differ- 
ent aspects of localization of different types of waves, 
e.g., electron wave, sound wave, electromagnetic wave, 
etc j and quantum matter wave in the form of a 
Bose-Einstein condensate (BEC) In the case of 

a quantum matter wave a disordered speckle potential 
and a quasi-periodic bi-chromatic optical lattice (OL) 
potential (| have been employed for the localization of 
a non- interacting cigar-shaped BEC. The quasi-periodic 
bi-chromatic OL potential was formed by the superpo- 
sition of two standing-wave polarized laser beams with 
incommensurate wavelengths. Each standing-wave po- 
larized laser beam leads to a periodic OL potential. The 
non-interacting BEC of 39 K atoms was created @ by 
tuning the inter-atomic scattering length to zero near a 
Feshbach resonance 

The localization of a cigar-shaped BEC in a one- 
dimensional (ID) bi-chromatic OL potential and related 
topics have been the subject matter of several theoretical 
flfjl - fl9l ] and experimental studies. (Localized BEC 
states in the form of gap soliton in a mono-chromatic 
OL potential have also been observed [2(| and studied 
[2l| . However, there are fundamental differences be- 
tween the two types of localized states. The Anderson 
localized states in bi-chromatic OL potential appear pre- 
dominantly in the linear Schrodiger equation and are 
destroyed above a small critical repulsive nonlincarity, 
whereas the gap solitons appear only in the presence of 
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a repulsive nonlinearity. The linear system does not sup- 
port a localized gap soliton. The gap solitons appear in 
the band gap of the spectrum of the periodic OL poten- 
tial. On the other hand, the potential supporting the 
Anderson states must be aperiodic in nature.) After the 
pioneering experiments [5, 6] on the localization of a ID 
cigar-shaped non-interacting BEC, a natural extension of 
this phenomenon would be to investigate localization in 
a weakly-interacting binary (two-component) BEC mix- 
ture. 

In this paper, with intensive numerical simulations of 
the Gross-Pitaevskii (GP) equation, we study different 
aspects of localization of a cigar-shaped binary BEC in 
a ID bi-chromatic quasi-periodic OL potential — similar 
to the one used in the experiment of Roati et al. [_§. 
Although, most of the studies of Anderson localization 
were confined to non-interacting systems, in the present 
study on the localization of a two-component BEC we 
shall mostly consider a weakly-interacting system under 
the action of both inter-species and intra-species inter- 
action. If the potential V(x) is symmetric around origin 
at x = 0, e.g. V(x) = V(—x), for the non-interacting 
system the density profiles are symmetric around x = 0: 
|u(a;)| 2 = \u(— x)\ 2 , where u(x) is the wave function of 
the matter wave. But for the interacting system with 
inter- and intra-species interaction, this symmetry can 
be broken and different types of symmetry-broken states 
emerge and we shall be specially interested in the study 
of the formation of localized cigar-shaped BEC compo- 
nents with broken symmetry. 

The two components of the BEC could be two differ- 
ent hyperfine states of the same atom [22j or two different 
atoms j23| and the two confining bi-chromatic OL poten- 
tial for the two components could be the same or different 
[H EH EII • The bi-chromatic OL potential could be taken 
as the linear combination of two sine terms, both with a 
minimum at the origin (x — 0), when the confined BEC 
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components will have a maximum at the origin. The Di- 
chromatic OL potential could also be taken as the linear 
combination of two cosine terms with a maximum at the 
origin, while the confined BEC components will have a 
minimum at the origin. We shall consider two distinct 
situations of localized binary BEC: (i) two components 
under the action of the identical sine bi-chromatic OL 
potential while both densities have maxima at x = 0, 
and (ii) component 1 under the action of the sine bi- 
chromatic OL potential with a maximum of density at 
x = and component 2 under the action of the cosine 
bi-chromatic OL potential with a minimum of density at 
x = 0. In case (i) with weak inter-species and intra- 
species repulsions, the localized BEC components are 
overlapping with a maximum at x = 0. But under the 
action of a inter-species repulsion above a critical value, 
the two repealing BEC components separate and move 
to opposite sides of x = 0, thus breaking the symmetry 
around x = 0. It is this separation or splitting of the 
BEC components which we study. (Similar separation of 
gap solitons [2lj in a binary BEC mixture confined by a 
monochromatic lattice has been studied [U,[2^].) In case 
(ii) a different type of symmetry breaking emerges. In 
this case with weak intra-species and inter-species repul- 
sion or attraction, the localized BEC components have 
densities symmetric around x = 0. But under the action 
of an inter-species attraction above a critical value, the 
two attracting BEC components try to stick together as 
a pair of bright solitons [26| and the minimum- energy 
stable configuration is the one where both components 
move predominantly to the same side of x = 0, thus 
breaking the symmetry around x = 0. We study this 
symmetry breaking in details. In our study we use direct 
numerical simulation of the underlying GP equation us- 
ing the split-step Fourier spectral method. We also use 
a time-dependent variational analysis for an analytical 
understanding of the numerical results in case (i) above 
when both components have a maximum at x = with 
densities having quasi- Gaussian shapes. In this case cer- 
tain aspects of the splitting of the two components are 
also predicted correctly from the variational analysis. 

There have already been a number of theoretical and 
experimental studies on Anderson localization in differ- 
ent problems. The effect of a repulsive non-linearity on 
localization of light waves in photonic crystals was stud- 
ied experimentally [27j ■ There have also been studies of a 
weak repulsion on localization [l2|, [TJ, [Tt], HH . Damski et 
al. and Schulte et al. considered Anderson localization in 
disordered OL potential [ll[ whereas Sanchez-Palencia et 
al. and Clement et al. considered Anderson localization 
in a random potential (Toj . There have been studies of 
Anderson localization with other types of disorder [29| . 
Anderson localization in BEC under the action of a dis- 
ordered potential in two and three dimensions has also 
been investigated (30j . 

In Sec. fll] we present a brief account of (i) the two- 
component ID GP equation, (ii) the bi-chromatic OL 
potentials used in our study and (iii) a time-dependent 



variational analysis of the GP equation under appropri- 
ate conditions. In Sec. IIIII we present our numerical 
studies on localization using the split-step Fourier spec- 
tral method. The density profile of the quasi-Gaussian 
localized states are in agreement with the variational re- 
sults. The density of the localized states are symmetric 
around the origin at x = in the absence of inter-species 
interaction. We study the symmetry breaking under the 
action of an inter-species interaction. Certain aspects of 
symmetry breaking can be explained with the variational 
analysis. We demonstrate the dynamical stability of the 
two-component localized states. In Sec. IIVI we present a 
brief discussion and concluding remarks. 



II. ANALYTICAL CONSIDERATION OF 
LOCALIZATION 

We consider a two-component cigar-shaped BEC un- 
der tight transverse harmonic confinement with the bi- 
chromatic OL potential acting along the axial x direction. 
Then it is appropriate to consider a ID reduction [3l[ of 
the 3D GP equation by freezing the transverse dynamics 
to the respective ground state and integrating over the 
transverse variables. Such a binary BEC is described by 
the following coupled ID BEC equation HE2: 

du\ 1 d 2 u\ 2 2 

i-gj- = ~2~Qx T +9l \ Ul \ +S' 12 M Ul + Vi[x)ui, (1) 

duo 1 2 Un , l9 , ,9 < . . . 

i^- = --^r4 +£2|u 2 | 2 +gi 2 K| 2 w 2 + V 2 (x)u 2 , (2) 



dt 



2 dx 2 



with normalization 



2 dx = 1, of the localized 



wave function uj of the two components j = 1,2. As 
the localization of the BEC in the bi-chromatic lattice 
is most prominent for the linear problem, and we shall 
be interested in small non-linearities gi,g 2 and gi 2 . To 
make the parameters of the model tractable we take the 
number of atoms Nj and the mass nij of the two com- 
ponents to be equal: Ni = iV 2 = N and m\ = ?n 2 = m. 
This simplification will have no effect on the general con- 
clusion of this study. The intra-species nonlinearities are 
given by (UJ gj = 2ajN/aj_,aj_ = \J % / (muj) , uj = y/ui y ui z 
where uj y and lo z are the transverse trap frequencies in 
the y and z directions, respectively, and Oj is the intra- 
species atomic scatte ring length. The inter-species non- 



linearity is given by }31| <?i 2 = 2Nai 2 /a±, where et 12 is 
the inter-species scattering length. In Eqs. Q and © 
Vj(x) are the bi-chromatic OL traps of the two compo- 
nents j — 1,2. Here we are using harmonic oscillator 
units: length is expressed in units of transverse oscillator 
length a± and time in units of angular frequency 

In the first part of our study we shall take the two 
bi-chromatic OL potentials to be identical — V\ (x) = 
V2(x) = V(x) — and have the following sine form as in 
the experiment Q 



V(x) = ^^sin 2 (fc ; a;), 



(3) 



i=i 
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with Ai — 2ir 2 si /Xf , where A; 's are the wavelengths of the 
laser forming the OL potentials, sj's are their intensities, 
and ki = 2it/\i the corresponding wave number. In the 
second part of our study we shall take V\ (x) to be given 
by Eq. ([3]) whereas in V 2 (x) the sine term is replaced by 
a cosine term: 



[321 ]. We shall obtain useful relations among the varia- 
tional parameters of the localized BEC states. This will 
allow us to get a better physical understanding of the 
localized states. 

The GP equations ([I]) and ([2]) can be derived from the 
following Lagrangian 



V 2 (x) 



i=i 



Ai COS 2 (fc;x), 



(4) L(t) = 



Potentials (j3j and ((4]) are quite similar. However, po- 
tential ((4]) generates a different type of localized states 
compared to potential ©. Potential ^ has a local min- 
imum at the center x = 0; consequently, stable stationary 
solutions with this potential have a maximum at x = 0. 
However, potential ((4]) has a local maximum at x = 
corresponding to a minimum of the stationary solution 
at the center. 

With a single periodic potential of the form 
sin 2 (27rx/A) or cos 2 (27rx/A) with s 2 = 0, the linear 
Schrodingcr equation permits only de-localized states in 
the form of Bloch waves. Localization is possible in the 
linear Schrodinger equation due to the "disorder" in- 
troduced through a second component in Eqs. © or 
(|4]). The primary lattice is usually strong enough and 
is used as the main periodic potential fixing the Bloch 
band structure of the single-particle states without dis- 
order. The secondary lattice is weaker (s 2 << si) and 
introduces a "deterministic" disorder. Following the ex- 
periment of Roati et al. Q, the transverse harmonic- 
oscillator length is of the order of a few microns so that 
the wavelengths of the OL in dimensionless units become 
approximately 10° ~ 10 1 . 

Usually, the stationary BEC localized states formed 
with bi-chromatic lattices occupy many sites of the quasi- 
periodic OL potential and their shape is modulated by 
the short-wavclcngth potential 



V(x)\u 
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gi 2 \ui\ 2 \u 2 \" 



dx, 



(6) 



where the overhead dot denotes time derivative, the star 
denotes complex conjugation, and the prime denotes 
space derivative. Using ansatz ([5]) in Eq. © we get 
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gi 2 \ui\ 2 \u 2 \ 2 dx 
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gi 2 NiN 2 



y/n{w 2 + W 2 ) 



exp 



Wn 



L ij = ^7 [cos(2/c ; x 0j )exp(-fc 2 w 2 ) - l] 



(7) 



(8) 



(9) 



where p = xq 2 — xoi , and the variational parameters are 
the norm Afj , width Wj , position xoj , and phase param- 
eters Cj and cj>j . We use the Euler-Lagrangc equation 

dL d dL 
da dt da 



0. 



(10) 



the parameters, potential leads to bound states con- 
fined practically to the central site of the quasi-periodic 
OL potential. When this happens, a variational approx- 
imation with Gaussian ansatz leads to a reasonable pre- 
diction for the bound state. To apply the variational ap- 
proach with potential §3§ effective on both components, 
we ado pt t he Gaussian function below as the variational 
ansatz 1331 



32j. For certain values of where a is the variational parameter. The first vari- 



ational equation using a 



yields Afj = constant. 



Without losing generality we shall take this constant to 
be unity and use it in the subsequent equations. Tak- 
ing a = xoj , Cj , Wj , and Afj , respectively, we obtain the 
following equations 



d 



Uj{x,t) 




^exp[-(x-x 0j ) 2 /(2^)] 



x exp[j{Cj(x - x j) + (j)j}\, 



dx j 
xoj = Cj , 
1 ^ Jh 



E Ll i + Ll2 



i=i 



(11) 
(12) 



(5) 



2ffi 2 w? 



where Wj are the widths of component j of localized BEC 
centered at xoj, <fij is the corresponding phase, Cj is 
the linear phase coefficient and Afj is the normalization 
of component j. We shall consider solution where the 
widths are comparable to the OL wavelength. Generally, 
the variational approach is valid when the OL is only 
smooth and slowly varying on the localized states scale 



V2^ 
x exp 



y/n(wl + w 2 ) 3 / 2 



2p 2 



^2A lW pf cos{2kix 0j ) 



1=1 



x exp(-kfw 2 ) = 0, 



9j 



Aw 2 



Wj V 27T 



dAfj 



Mj=l 



(13) 

E^- 

i=i 

(14) 
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The set of equations (JTT|) — (fl4j) shows that the phase 
<pj has no effect on the location xqj, width Wj and the 
coefficient Cj of solitons, whereas it is determined by 
the variational equation (|T4"|). Hence we shall neglect 
Eq. (jT4j) in the following and consider only Eqs. — 
(fT3|) . (However, we could not set phase <f>j — in the 
beginning, as the Eulcr-Lagrange equation for <pj leads 
to the important result of norm conservation.) Equation 
(fT3l) determines the width wj in terms of nonlincaritics 



9j: 9 12 an d the parameters of the bi-chromatic OL po- 
tential. The density distribution is quasi-Gaussian only 
for small values of xoj, and we shall be limited to this 
constraint (small \xoj\) while considering the variational 
approach. In the symmetric case, while gi = g2, the 
widths of the two localized states are equal. However, 
when c/i 7-52, the widths of the two localized states are 
not equal and one has two asymmetric localized states. 
When xqi = X02 = 0, the widths of the two stationary 
BEC localized states are determined by 



fJj 



Zgiiw? 



Wj 



2ir V^( w l + wf) z/2 



1=1 



2Aiw^kfexp(-kfw 2 ) = 



(15) 



Now we obtain, following Rcf. [34(, the equation of 
motion describing the center of the two BEC localized 
states as if they were a particle in an effective potential. 
Thus, we can determine the motion of equivalent parti- 
cles using the classical mechanics analogy. Inserting Eqs. 
©, © and (TTTl) into Eq. (TT2"]l. we obtain the following 
equation for dynamics 



d xoj 

v cS3 



dV c sj 
dx j 
gi2 



(16) 



2 A, 



■. exp 



- E y[cos(2fc,x,y) exp(-fc^ 2 ) - 1] , (17) 
1=1 

where the an-harmonic effective potentials V B sj are the 
dynamic potentials for the movement of the center of the 
localized BECs and depend on characteristics of the two 
BEC localized states (e.g., width wj and position xoj) 
and potential parameters (e.g., ki and A{). The effective 
potentials have two parts. The first term on the right of 
Eq. (|T7|) is induced by the mutual interaction of the two 
BEC localized states. When the coupling constant g\2 > 
0, it is positive and contributes to a repulsive potential 
barrier for BEC localized states. The second term in 
Eq. (|17|) arises from the bi-chromatic OL potential and 
contributes to an attractive potential well, if \xoj \ is small 
enough. The combination of the two terms may lead to 
local minima or local maxima at the origin xoj = in 
the effective potential V e sj- 

When the two BEC localized states are symmetri- 
cal (gi = 172), they will feel the same effective poten- 
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FIG. 1: (Color online) (a) The effective potential 14ff(a;o) vs. 
the position of the localized state xq = xoi = — £02 in the 
symmetric case gi = = for 1712 = 0.5, 1 and 3 as obtained 
from Eq. (|17|) . The respective widths are calculated from Eq. 
(|13p . (b) The effective potential felt by the first component 
Veffi(xoi) vs. its position xoi as obtained from Eq. (|17|) . 
The widths of the two states are calculated from Eq. (|15|) . 
In both cases — (a) and (b) — the potential parameters are 
Ai = 8, A 2 = 0.862Ai, si = 8, s 2 = 2.4. 



tial. In that case, we should have w± — W2 — w and 
xqi = —X02 = xq. The effective potential V c s is plotted 
in Fig. [T] (a) vs. xq in the symmetric case g\ = .92 = 
for different inter-species interaction: g\2 — 0.5,1, and 
3. Equation (|13[) is first solved to find w as a function 
of xq. This result is then substituted in Eq. (fTT)) to 
find V4ff(a;o) as plotted in Fig. [T] (a). From Fig. Q] (a) 
we find that if g\2 is smaller, the effective potential pos- 
sesses a local minimum at xq = 0. In this case if we start 
with two BEC localized states at x = for small 312, 
they will be in stable equilibrium and unsplit. The effec- 
tive potential becomes weaker with the local minimum 
at xq — less pronounced as gi2 increases. If 1712 is large 
enough (e.g., 1712 = 3 in Fig. [T] (a)), a local maximum 
appears in the effective potential at xq = 0; and the peak 
value at maximum increases as g±2 increases. In this case 
two overlapping (unsplit) initial localized BEC states at 
the maximum at x = will be in a metastable configu- 
ration. When they arc slightly perturbed, for example, 
by slightly moving the BEC center(s) from x — 0, they 
will move respectively toward the minima of the effective 
potential next to the maximum and split into two non- 
overlapping BECs localized at the two local minima of 
the effective potential on both sides of xq = as we shall 
see in Sec. IHI AI 

When the two BEC localized states are asymmetrical 
(.9i 7^ .92), they will feel different effective potentials, ft 
is difficult to calculate these effective potentials exactly. 
But we can calculate the effective potential(s) (fT7| un- 
der simplifying assumptions. The widths w\ and W2 are 
assumed to be independent of positions xqi and xq2 and 
were calculated using Eq. (fTS"]) and substituted into Eq. 
(fTT)) . The function V r ffi(xoi) so obtained is plotted in 
Fig. (TTJ) (b) vs. x 01 , for g x = 0,g 12 = l,g 2 = and 2. 
The effective potential acting on component 1 is stronger 
as 172 is increased from to 2. This will lead to a reduc- 
tion of the width of the corresponding localized state as 
172 is increased from to 2 as we shall see in Sec. III! Al 
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III. NUMERICAL RESULTS 

We perform the numerical simulation employing the 
real-time split-step Fourier spectral method with space 
step 0.04, time step 0.001. (We also checked the re- 
sults using the real-time Crank-Nicolson discretization 
routine, the FORTRAN programs which are given in Ref. 
(3f|.) Because of the oscillating nature of the potential, 
great care was needed to obtain a precise localized state. 
The accuracy of the numerical simulation was tested by 
varying the space and time steps as well as the total num- 
ber of space steps. The initial input pulse is a Gaussian 
wave packet, Uj(x,t = 0) = 7T" 1 / 4 exp(— x 2 /2) with a 
parabolic trap V'(x) — x 2 /2. Let g\ — gi = 912 = 
and the coupled Eqs. ((T|) and © become the linear 
Schrodingcr equation. To start the numerical simula- 
tion the parabolic trap is slowly turned off and the bi- 
chromatic OL is slowly turned on; the increment in the 
coefficient s, in each time step is 0.00005. Successively, 
we add gradually the nonlinear coefficient 32 and the cou- 
pling constant 312 slowly. The increment of gi and gn is 
0.00001 in each time step. Thereby the stationary local- 
ized states are obtained. 

To perform a systematic numerical study of localiza- 
tion we take s\ = 8, si = 2.4, Ai = 8 and A 2 = 0.862Ai 
maintaining the ratio A2/A1 = 0.862 similar to the value 
employed in the experiment 0. The values of si and S2 
are also in the range employed in the experiment. With 
these set of parameters and for weak inter-species and 
intra-species interactions, the localized BEC state could 
be mostly confined in the single cell of the bi-chromatic 
OL potential. All results reported here are obtained with 
these potential parameters except those in Figs. [3j(e) and 
(f) where we study the effect of the change of potential 
parameters. We consider a transverse harmonic oscilla- 
tor length a± ss 1 /im, so that the wave lengths of the 
OL potential arc approximately 800 nm and 689 nm. We 
shall take g\ = throughout the present investigation. 

A. Localized states with potential ([3]) on both 
components 

We start the numerical analysis with a consideration of 
the widths of the two components of the localized BECs 
with potential (|3J) acting on both components. In this 
case both the components could have a maximum at the 
center x = with localized states of quasi- Gaussian forms 
and when this happens the system can be described well 
by the variational approximation. Assuming a Gaussian 
distribution for the localized BEC, the numerical width 
is calculated via 

/oo 
\ Uj \ 2 x 2 dx. (18) 
-00 

The numerical and variational results for the widths Wj 
vs. <?2 for gi2 = 1 are plotted in Fig. [5] (a). The same 
for 172 = 0.5 vs. <7i2 are plotted in Fig. H] (b). In these 




(a) (b) 

FIG. 2: (Color online) (a) The numerical and variational 
widths wi and u>2 of the two stationary BEC localized states 
vs. Q2 for (?i2 = l,<?i = 0. (b) The same as in (a) vs. 1712 
for Q2 = 0.5, (71 = 0. The potential (UJ acts on the two com- 
ponents where the parameters are Ai = 8, A2 = 0.862Ai , si = 
8,s 2 = 2.4. 



plots there is reasonable agreement between the numer- 
ical and variational results for the width and the width 
usually increases with the increase of nonlinearities gi 
and 312 except in the case of w\ in Fig. [5] (a). The gen- 
eral increase of width with the increase of non-linearity 
(increase of repulsion) is expected, however, the decrease 
of w\ with gi in Fig. [2] (a) requires special attention. 
This phenomenon will be clear if we consider the effec- 
tive potential of component 1 as illustrated in Fig. [1] (b) 
when 32 is changed from to 2. We find that the effec- 
tive potential becomes slightly stronger as gi increases 
corresponding to a reduction of width w\ in Fig. [2] (a) 
with the increase of gi. 

Next we consider typical numerical profiles of the 
density of the localized BEC states |uj| 2 and compare 
them with variational results for small values of 312 
when the centers of both components are at the origin: 
a^oi = £02 = 0. In Figs. [3J (a) and (b) we present 
the results for g\ = #2 = and g\2 = and 1, re- 
spectively, while the two density profiles |ui| 2 and |i*2 1 2 
are equal. Next we consider results for the asymmetric 
case while g\ 7^ 92 while the two densities are differ- 
ent. In Figs. [3j (c) and (d) we present the results for 
.91 = 0,52 = 1,512 = 2 and for gi = 0,.g 2 = 2,312 = 1, 
respectively. In both cases the component 1 is more 
strongly localized in space with a smaller width. Next 
we consider the effect of a variation of the wavelengths 
Ai and A 2 on the densities. In Fig. [3J (e) we plot the 
densities for Ai = 10, A 2 = 0.862Ai, and in Fig. [3J (f) we 
plot the densities for Ai = 5, A2 = 0.862Ai; in both cases 
we take g\ = 0,gi = .912 = l> s i = 8, s 2 = 2.4. In Fig. 
[3J (e) with the increase of A's the central OL site has a 
larger spatial extension, consequently the densities have 
a smaller value at the maxima at the origin correspond- 
ing to a larger width. In Fig. [3J (f) with the decrease of 
A's the central OL site has a smaller spatial extension, 
consequently the densities have a larger value at the max- 
ima at the origin corresponding to a smaller width. The 
plots in Figs. [3J (a) — (e) indicate that the numerical 
results for densities are in perfect agreement with the 
variational results. It also indicates the good accuracy 
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FIG. 3: (Color online) The densities |wi(a;)| 2 and |u2(a:)| 2 
of the two localized BEC states vs. x for (a) 3i = 32 = 
312 = 0; (b) 31 = g 2 = and 312 = 1; (c) 31 = 0,32 = 1, 
and 312 = 2; (d) 31 = 0, 32 = 2, and 312 = 1. In (a), (b), 
(c), and (d) potential (O acts on the two components with 
parameters Ai = 8, A2 = 0.862Ai,si = 8, S2 = 2.4. The same 
as in (a) for 31 = 0, 32 = 312 = 1 with potential parameters 
(e) Ai = 10, A 2 = 0.862Ai,si = 8, s 2 = 2.4 and (f) Ai = 
5, A2 = 0.862Ai,si = 8, S2 = 2.4. In each case the profile of 
the bi-chromatic OL potential is also shown. 



(a) (b) 

FIG. 4: (Color online) (a) The same as in Fig. [3] (a) for (a) 
31 = 32 = and 312 = 3, and (b) 31 = 0, 32 = 0.3 and 
312 = 3. The parameters of the potential are Ai = 8, A2 = 
0.862Ai,si = 8, s 2 = 2.4. 
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FIG. 5: (Color online) The phase diagram of 312 vs. 32 for 
31 = showing the split and unsplit configurations of the 
two localized states while potential ((3| is effective on both 
components. The parameters of the potential are Ai = 8, A2 = 
0.862Ai,si = 8, s 2 = 2.4. 



of the numerical routine. In these cases the density has 
a quasi-Gaussian profile used in the variational analy- 
sis. In Fig. [3] (f) the numerical densities have secondary 
maxima on both sides of the central peak, and the profile 
deviates from the Gaussian and the agreement with the 
variational results in this case is only fair. 

So far we presented results in the asymmetric case for 
small values of 312 where the two localized states are over- 
lapping (unsplit) and centered at x = 0. As 312 increases 
it is clear from Fig. [T] (a) that the effective potentials 
experienced by the components develop a maximum at 
the origin. Hence the position x = ceases to be one 
of stable equilibrium for the components. Consequently, 
the two localized BECs move on opposite sides of the 
point x = to attain split (separated) configurations of 
stable equilibrium. We illustrate this in the symmetric 
case in Fig. 2] (a) for gi = g 2 = 0,512 = 3 where we plot 
the densities of the two split components. The two densi- 
ties continue to be symmetric: |wi(x)| 2 = |u2(— x)| 2 with 
centers at |xoi| = |xo2 1 = 0.45. From Fig. [T] (a) we find 
that the minima of the effective potential for the same set 



of parameters as in Fig. 2] (a) are at |xoi| = |xo2| = 0.58 
indicating the centers of the split solitons to be at ±0.58. 
The variational analysis is valid for small splitting, and 
considering that the value ±0.58 indicating splitting is 
not too small compared to the extension of the localized 
states, the agreement between the numerical displace- 
ment 0.45 of the localized state and its variational esti- 
mate of 0.58 should be considered satisfactory. In Fig. [4] 
(b), we present an example of the localized asymmetric 
split states for gi = 0,32 = 0.3, 312 = 3, where the two 
states have different (asymmetric) densities. 

As the split and the unsplit configurations of the local- 
ized states are of interest, it is appropriate to represent 
the split and unsplit configurations in a 32 vs. g\ 2 phase 
diagram for 31 = 0. This is illustrated in Fig. [5j (Sim- 
ilar split-unsplit phase diagram has also been studied in 
the context of symbiotic gap solitons 24[.) For small 
312, as expected, the localized states are in an overlap- 
ping (unsplit) configuration. The split configuration is 
achieved when 312 is larger than a critical value indicat- 
ing a minimum of intcr-specics repulsion needed to move 
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FIG. 6: (Color online) The densities |mi(:e)| 2 and |«2(x)| 2 of 
the two localized BEC states vs. x for (a) gi — g 2 = g\ 2 = 0, 
(b) gi = g 2 = 0,312 = -4, (c) .gi = 0, g 2 = -0.5, g 12 = -3, 
and (d) gi = 0, g 2 = 0.5, gi 2 = —3, when potential ^ acts 
on the first component and potential @ acts on the sec- 
ond component. The potential parameters are Ai = 8, A2 = 
0.862Ai,si = 8, s 2 = 2.4. 
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FIG. 7: (Color online) The phase diagram of \gi 2 \ vs. g 2 for 
gi = showing the symmetric and asymmetric configurations 
of the two localized BEC states. The potential © is effec- 
tive on the first component and potential on the second 
component with the same parameters as in Fig. 



the localized states from x = 0. Splitting is favored when 
the inter-species interaction is strongly repulsive with a 
large positive value of g\2 and when the intra-species in- 
teraction is attractive corresponding to a negative g2 as 
seen in Fig. [5] 



B. Localized states with potential ([3]) on 
component 1 and potential (|4]) on component 2 

Now we consider the situation when potential © acts 
on component 1 and @ acts on component 2. In this 
case only the first component under the influence of po- 
tential @ could have a maximum at the origin (x = 0), 
and second component under the influence of potential 
(U) usually has a minimum at the center and hence the 
variational approximation is not applicable in this case. 
(The monochromatic counterpart to this model with the 
phase shifted potential (@| was investigated in Ref. [3q|.) 
We start our discussion in this case by showing results of 
density under different situations. In Fig. [6] (a) we plot 
the densities of the two components for the symmetric 
non-interacting case for g\ = g^ = g\2 = 0. The density 
I ui 1 2 of the first component has a localized Gaussian- 
type shape centered at x — 0. The density of the second 
component has a two-hump symmetric structure with 
a minimum at x = 0. The symmetry of the densities 
around x = can be broken in the presence of a suffi- 
ciently strong inter-species attraction. This is illustrated 
in Fig. |5] (b) where we plot the densities for gi = g 2 = 
and g\2 = —4. The center of both the densities has 
moved slightly towards positive x as a sufficiently strong 



inter-species attraction is introduced via .912 = —4. In 
Figs. [U (c) and (d) we illustrate the densities in two 
more cases for g\ =0,(72 = —0.5, and 512 = —3 and for 
gi = 0, 92 = 0.5, and g±2 = —3, respectively. In the case 
of Fig. |6](c), the densities are asymmetric and in the case 
of Fig. |5] (d), they are symmetric. This observation is 
consistent with the general conclusion that the symmetry 
breaking will be favored when the system is attractive. 
In the case of Fig. [5] (c) we have g± = 0,g2 = —0.5, 
and gi2 = —3 indicating that both components could 
be bound due to inter-species attraction even in the ab- 
sence of the bi-chromatic OL potential. Hence due to 
the inter-species attraction the two localized BEC states 
tend to overlap and stay together as two bright solitons 
[2r| . In the case represented in Fig. [5] (d) g 2 is positive 
and the system is less attractive and hence the density 
distribution is symmetric. 



Again it is worthwhile to show the symmetry breaking 
in a g2 vs. \gi2 \ phase diagram for g\ = as illustrated in 
Fig. [7] (Similar symmetry breaking has also been studied 
in the context of symbiotic gap solitons (24[.) For small 
I512I, as expected, the localized states are in a symmetric 
configuration. The asymmetric configuration is achieved 
when I <7i 2 1 is larger than a critical value indicating that 
a minimum of attraction is needed to replace the center 
of the localized states from x = 0. Symmetry breaking is 
favored when the inter- and intra-species interactions are 
strongly attractive leading to strong effective attraction 
in the system and hence favored for large negative (at- 
tractive) values of 172 than for positive (repulsive) values 
of g 2 as seen in Fig. [71 
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(a) (b) 

FIG. 8: (Color online) The densities (a) \u!(x,t)\ 2 and (b) 
|lt2(as, t)| vs. x and t of the two localized BEC states of Fig. 
[4] (b) when the potential parameters are suddenly changed 
from Ai = 8, A 2 = 0.862Ai, si = 8, s 2 = 2.4 to Ai = 8.5, A 2 = 
0.862Ai, si = 8, s 2 = 2.4 at time t = 5. 



C. The stability of the localized BEC states 

So far wc studied the stationary properties of the local- 
ized BEC states. Now we investigate if the localized BEC 
states are stable under small perturbation. We are using 
real-time propagation routines in the numerical calcula- 
tion, which finds usually the stable states. This indicates 
that the localized BEC states, we are studying, should be 
stable. However, we now explicitly demonstrate the sta- 
bility of a set of specific states, e.g., the ones illustrated in 
Fig. |4](b). To demonstrate that these states are stable, 
after their formation, we suddenly change the potential 
parameters from Ai = 8,A2 = 0.862Ai, si — 8,s-2 = 2.4 
to Ai = 8.5, A 2 = 0.862Ai, si = 8, s 2 = 2.4 and study the 
resulting dynamics. We plot in Figs. [8] (a) and (b) the 
consequent density dynamics |ui(:c,£)| 2 and |M2(a;,i)| 2 vs. 
x and t, respectively, for t = to 20. The sudden change 
of the potential is effected at t = 5. From Figs. [5] 
(a) and (b) we see that at t = 5 when the potential 
is suddenly changed the density profiles suffer an abrupt 
change. Nevertheless, the BEC density profiles remain 
localized for a large interval of time undergoing breath- 
ing oscillation as we can see from Figs. [8] (a) and (b) for 
the two components, which demonstrates the stability of 
the localized BEC states. 



IV. SUMMARY 



a maximum of density at the trap center (a; = 0) satisfy- 
ing the symmetry |uj(a;)| 2 = \uj(—x)\ 2 . In this case the 
variational analysis is applicable and produced result in 
good agreement with numerical simulation. In case (ii), 
for a weak inter-species interaction, one component has 
a maximum of density at the trap center and the other 
component has a minimum, however, again satisfying the 
symmetry |uj(a;)| 2 = \uj{— x)\ 2 . The localization is most 
favored for non-interacting atoms and is destroyed in the 
presence of moderate intra-species and inter-species re- 
pulsive interactions. Here we specially study the effect of 
weak inter-species and intra-species nonlinearities on the 
profile of the localized states. In case (i), for a sufficiently 
strong inter-species repulsion, the two localized states are 
found to move in opposite directions from the trap center 
and attain equilibrium in a separated (split) configura- 
tion breaking the symmetry |uj(a;)| 2 = \uj(— x)\ 2 . Incase 
(ii), for a sufficiently strong inter-species attraction, the 
two localized states are found to move away from x = to 
the same side of the trap center again breaking the sym- 
metry \uj(x)\ 2 = \uj(— x)\ 2 . Wc studied in some detail 
the formation of the localized states of broken symme- 
try in all cases considering phase diagram of intra- and 
inter-species nonlinearities gi and gn- 

The Anderson localization of a BEC in a bi-chromatic 
OL potential as studied here is very distinct from the 
localization of gap solitons [2(| [U l36j of the BEC in 
a mono-chromatic OL potential. Anderson localization 
takes place in an aperiodic potential in the linear equa- 
tion, whereas a periodic potential and a repulsive nonlin- 
earity are crucial for the localization of gap solitons. 

We hope that the present work will motivate new stud- 
ies, specially experimental ones, on the localization of a 
binary mixture of BEC in a bi-chromatic OL potential. 
The effect of nonlinear inter- and intra-species interac- 
tions on localization, as investigated in the present paper, 
also deserves careful analysis. 
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In this paper, using the numerical and variational so- 
lution of the time-dependent CP equation, wc studied 
the localization of two-component cigar-shaped BEC in a 
quasi-periodic ID OL potential prepared by two overlap- 
ping polarized standing-wave laser beams with different 
wavelengths and amplitudes. Specifically, we considered 
two analytical forms (sine and cosine) of the OL poten- 
tial and considered two distinct cases: (i) the two com- 
ponents under the action of the sine potential ([3]) and (ii) 
the sine potential §3§ acting on component I and the co- 
sine potential (0} acting on component 2. In case (i), for 
a weak inter-species interaction, both components have 
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